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Langevin/Fokker-Planck processes can be immersed in a larger frame by adding 
fictitious fermion variables. The (super) symmetry of this larger structure has been 
used to derive Morse theory in an elegant way. The original physical diffusive mo- 
tion is retained in the zero-fermion subspace. Here we study the subspaces with 
non-zero fermion number which yield deep information, as well as new computa- 
tional strategies, for barriers, reaction paths, and unstable states - even in non-zero 
temperature situations and when the barriers are of entropic or collective nature, as 
in the thermodynamic limit. The presentation is self-contained. 



Many systems have a dynamics with processes which take place on distinct timescales. 
The most familiar example is the diffusion in a many- valley energy E{x) landscape at low 
temperatures: 



which consists of rapid gradient descents into local minima, and slow 'activated' transitions 
between minima induced by the thermal noise rji (Gaussian independent white noises of 
unit variance). The relevant parameter is the inverse temperature (3 = 1/T: the larger 
f3 the more pronounced the gap between fast intra- valley relaxations and slow activation 
processes ln{tsiow) ~ P^E Another example is that of cooperative systems at finite 
temperature. Consider for example a ci-dimensional ferromagnet in a magnetic field h point- 
ing upward: the state with negative magnetization becomes unstable, its decay taking a 
time \ia{tsiow) oc h~^'^~^\ the parameter controlling the timescale separation is the inverse 
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of the field. In the absence of field, the slow relaxation takes a time ln(ts/oto) oc L'^~^, and 
the control parameter is the size L}^. From a conceptual point of view, it is important in 
these situations to characterize the relevant structures: metastable states and their basins 
of attraction, the reaction paths joining them, and the timescales involved. On the other 
hand, in order to efficiently model realistic situations, one needs to be able to treat the rare 
'activated' passages in a specific way. 

Quite generally, a separation between timescales implies the existence of metastable 
states, defined as probability distributions corresponding to situations in which everything 
fast has happened and everything slow has not taken place. If there are more than two 
timescales, for example tjast ^ iinterm tsiow one has metastable states at tinterm, and a 
different set of metastable states at tsiow, the latter resulting from the fusion or the decay of 
states defined for tinterm- Given how natural the concept of metastability is, it may come as 
a surprise that only recently has a construction of metastable states based on the stochastic 
dynamics been fully estabhshed j^. 

The idea is simple: the probability associated with (|TJ evolves through jl^ : 



(from here onward we denote derivatives as A^i = ^ and A^ij = gf^r) where Hpp is 
the Fokker-Planck operator. It turns out that if there is a separation of timescales 
tfast tsiow in the system, the spectrum of Hpp has a gap: there are (say) K eigenvalues of 
the order of tjj^^ and all other eigenvalues are at least of the order of tj^^^. Furthermore, one 
can show that, in the limit of large timescale separation, irrespective of its origin, one can 
construct exactly K probability distributions corresponding to distinct metastable states by 
linear combinations of the K eigenstates 'below the gap'. The low temperature example is 
particularly clear: in that case these K distributions are Gaussians of width ^/T sitting at 
each of the K local minima. 

In a situation with metastability, it becomes interesting and in practice necessary to 
evaluate the time of decay, as well as the spatial probability distribution of the escape 
current (i.e. the reaction paths). This may involve identifying the barrier or 'bottleneck' 
responsible for the slowness of decay. In the low temperature example, the reaction paths 
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are simply gradient lines connecting two energy minima through a saddle with one unstable 
direction (of index one). The bottlenecks are these saddle-points , an d there is a considerable 
variety of methods for their location in high dimensional space (l5l |. 

Now, there is a construction that naturally incorporates saddle points, and that has been 
successfully used to derive the relations between the numbers of saddle points of a function 
(the energy here) and the topological properties of the manifold on which it is defined. These 
relations are the so-called Morse inequalities , and have been rederived in an elegant and 
elementary way by Witten The construction is a 'completion' of the diffusive problem 
as follows. First express the Fokker-Planck operator in a basis in which it is manifestly 
Hermitian jl^ : 
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Second, 'complicate' the operator and the space by introducing N fermions a|, Cj, with 
z = 1, . . . , A^, and: 



Within zero-fermion subspace, is just the Hermitian form of the Fokker-Planck operator, 
and, as we have remarked, for low temperatures its eigenstates 'below the gap' are related 
to local minima. 

The wavefunctions with one or more fermions are so far a spurious addition. Their interest 
stems from the fact that one can show that there is a gap in all the spectra associated 
with any fermion number, and the states 'below the gap' having one fermion are for low 
temperatures peaked on saddles with one unstable direction, those having two fermions on 
saddles with two unstable directions, and in general those with p fermions on saddles of index 
p. Using this fact and the symmetries of H^, it is then easy to derive Morse inequalities 
- we shall review this in Section 11111 

The basis that makes Hpp Hermitian as in (jS)) (we shall in what follow refer to it as the 
'Hermitian basis') offers the direct way to get to Morse theory, and is the one most often used 
in the field theory literature. Here, we are mainly interested in the diffusive interpretation, 
at least of the original zero-fermion subspace, and for this we must go back to the original 
basis The change from the Hermitian to the original basis is made by a multiplication 
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by e'^^/^, and is quite tricky since it is exponentially large in the relevant parameter — /3 if 
we are interested in the low temperature limit, or the system size for macroscopic systems. 
Because of this reason, one has to be very careful because negligible, large deviations in one 
basis become of 0(1) in the other. Indeed, the strategy we shall follow in this paper is to 
rederive the limit wavefunctions in each basis from scratch. 

A first question we may ask is how do the eigenstates 'below the gap' with fermion number 
larger than zero look, for low temperatures, in the original basis in which 

H = Hfp + J2 E,,a]a, = e-^^^'H^e^^^^ (5) 

The outcome, as we shall see, is a pleasant surprise: for example one fermion (right) eigen- 
states 'below the gap' are concentrated not on the saddle, but along a narrow (width ~ VT) 
tube following the gradient line joining minima and passing through the saddle - the reac- 
tion path. Higher fermion number subspaces (and left eigenstates) also encode interesting 
information. 

The construction yielding Morse theory relies on the low-temperature limit, in which 
functions peak on the appropriate structures. Low temperatures are just one instance of 
evolution with widely separated scales. One is naturally led to ask what happens with 
the construction we have described in the presence of a timescale separation generated by 
some other (collective, entropic.) mechanism. Again, the answer is pleasant: for example 
the one fermion eigenstates 'below the gap' of ^ yield the reaction current distributions 
between metastable states (the latter defined dynamically as outlined after Eq. Q). This 
generalization will give us practical strategies for the evaluation of reaction paths, valid 
whenever there is timescale separation. From a more abstract point of view, it will yield 
a precise definition of 'free energy barrier' in a natural way, without having to rely on 
mean-field or any other approximation. 

Let us write, for a generic wavefunction Itf)), an evolution equation: 

dm 



dt 



-H IV). (6) 



Specializing \j) to zero fermions we recover the Fokker-Planck equation Consider 
now (jH)) but within the one-fermion subspace, in which functions are of the form |^^) = 
j d^K Xlc -^c(^c)|x)aj| — ), with \x) the basis in space and | — ) the fermion vacuum. It 
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amounts to an evolution for a fecior function R{x,t) = {Ri{x,t), .., RN{x,t)): 



-^^-H,,KX.,t)-Y^—R,A.,i). (7) 

b 

Equation ^ is one of the main instruments of this paper. It evolves a vector field R{x^ t) 
so that it rapidly becomes a linear combination of one fermion states 'below the gap'. 

In a system with metastable states, starting from an initial condition, the probability 
distribution P(a;, t) evolves rapidly to a quasi- stationary distribution corresponding to quasi- 
equilibrium within one or more metastable states. At longer timescales, P{x,t) will be 
gradually concentrated on new, more stable metastable states. If we choose the initial 
condition to be close to a state, we get a quick thermalization within such a state. What we 
have been discussing up to now suggests that equation does for reaction currents what the 
Fokker-Planck equation does for states: depending on the initial conditions, R{x,t) tends 
rapidly to a reaction current between metastable states, which, in the particular case of very 
low temperatures, is a single reaction path. As time passes, other new current distributions 
start contributing to R{x,t). 

The main point is that while escape from a state takes by assumption long times for the 
probability distribution P{x,t), convergence to a reaction current is for R{x,t) immediate 
(or more precisely, of the same order of the time it takes for P{x, t) to stabilize in the closest 
state). Now, a whole set of practical methods, such as simulated annealing and transition 
path sampling, can be seen as ways of implementing the Fokker-Planck equation — the 
former in a diffusive and the latter in a functional way. The same can be done with equation 
([7j): we shall give in this paper a diffusion equation which is to (|7j) what the Langevin process 
is to (j2I), and a path-sampling procedure based on (jZj) which has the peculiarity that the 
paths pile up on the barriers. 

To conclude this rather technical introduction, let us summarize what we do in this 
paper. We first (section |n|) introduce the supersymmetric construction in detail, stressing 
in particular the relation between the original Fokker-Planck basis and the Hermitian basis. 
As a first example, in section ITTTl we rederive Morse theory for the case of smooth potentials 
in a simply connected space. 

In section IIV| after briefly reviewing the djTiamic definition of metastable states, we 
analyze in detail the structure of the one fermion subspace, showing that it contains the 
reaction paths and also 'loops': reactions leading from a state to itself. Although for clarity 
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we always keep the low-temperature case in mind, the developments are valid whatever the 
origin of the timescale separation. 

In section |3 we summarize the structure in all fermion subspaces. For the low- 
temperature case we discuss the form of all eigenstates 'below the gap'. The derivation 
of these results can be made in an explicit way using a diffusive dynamics we introduce in 
section IVIl which plays for the higher fermion number subspaces the role that Langevin 
equation plays for the zero-fermion Fokker-Planck evolution. 

In section IVTTl we construct a path-sampling process to locate reaction paths, in which the 
trajectories are weighed with the usual Langevin action plus a Lyapunov exponent associated 
to each trajectory. 



II. SUPERSYMMETRIC QUANTUM MECHANICS AND FOKKER-PLANCK 

EQUATION 

We assume that E{x) is two times differentiable and bounded so that the Gibbs measure 
exists: 

d^xe-^^<oo. (8) 



The dynamics of the system is given by the Langevin equation (^. The probability dis- 
tribution will evolve according to the Fokker Planck equation which can be seen as a 
continuity equation for the current 1^ : 

Mx, t) = (t-^ + eA P{x, t). (9) 



We introduce at this point fermion creation and annihilation operators a| and Oj, with 
anticommutation relations [aj,a^]+ = 6ij; the fermion number operator is Nj = Xli'^l'^j- 
We denote states in the coordinate space as and ipix) = {x\iIj), using the Dirac bra-ket 
notation of Quantum Mechanics; the zero-fermion state is |— ), and we denote states in the 
product space (coordinate fermions) with boldface. We say that has n fermions if: 

iV;|V)=n|V). (10) 

We define the 'charges' 



7 



which satisfy: 

= = 0. (12) 

We can write the operator 

H = ^{Q + Qf = i[g, g]+ = Hfp + Y, E,i,a]ai, (13) 

where 

[H,Q] = [H,Q] = Q. (14) 

H commutes with the fermion number operator Nf, so that eigenstates are classified accord- 
ing to their fermion number. Within the zero-fermion space, H is the original Fokker-Planck 
operator. These relations are true only for the Fokker-Planck equation with the drift forces 
at least locally a gradient. Q and Q commute with H and transform states with an even 
number of fermions (bosonic states) into those with an odd number of fermions (fermionic 
states) and vice- versa, hence the name 'supersymmetry'. What we have done up to now 
can be seen as completing the square, and making symmetries underlying the Fokker-Planck 
equation with gradient forces explicit. 

One can now make a change of basis such that the charges become Hermitian conjugates 
of one another: 



= e^^l'Qe-^^l' = -i V(rA _ ^E.^a,, 

^-^ oxi 2 

i 

= e^^l-Qe-^^l- = -i ^^(r A + ^_eM = (Q'Y ■ (15) 



In the new basis we have the Hermitian equivalent of (|13|): 



rj-] L ^ 7 ^ J ~r rj-] \ 

fJ2 + 1^' - + E ^-«]«^ = (H'V- (16) 



Hpp has now the standard form of a Schrodinger operator (acting in imaginary time) with T 
playing the role of h. On the other hand is the standard Hamiltonian of Supersymmetric 



Quantum Mechanics 



24|. 



As the original operators Hpp and H are not hermitian we will have two different eigen- 
value equations, one for the right eigenstates and one for the left eigenstate ((V'"^|) 
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(17) 



while in the Hermitian basis there will be only one equation 



(18) 



The three states are related by 



(19) 



It is clear that H and H have the same spectrum. Furthermore, the relation: 



(20) 



implies that a left fc-fermion eigenstate of if is an — right eigenstate of the problem 
with the inverted potential —E. 

From ()13p and we see that H and have non- negative eigenvalues. By construction, 
there is at least one eigenstate corresponding to the eigenvalue A = (the smallest 

possible): 



In order for to be normalizable we need the convergence of the integral 

The left and right eigenstates ()21|) have zero fermions and, thus, they belong also to the 
spectrum of Hpp. Clearly, both Q and Q annihilate li/?"^). It is easy to show (see Appendix 
that this is necessary for any zero eigenvector, and indeed (j2H) are the only ones with 
this property if the space has no holes. 

In general applying Q to any eigenstate we get either a degenerate eigenstate with 

one less fermion or zero. Similarly, applying Q we get either a degenerate eigenstate with 
one more fermion or zero. Each non-zero energy eigenstate annihilated by Q, can be 

written as = Q\x^)j ^"^^ the same holds for Q. Indeed, from the eigenvalue equation 



1'!/'°"^) OC constant 



(21) 



(22) 
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FIG. 1: The pairing of the energy levels for a generic spectrum. Each eigenstate of positive 
energy has a supersymmetric partner. The number of fermions is written below each corresponding 
column. The only unpaired state is the zero-energy one. 

one can infer that 

Ix"") = ^OIV'^), (23) 

satisfies Q\x^) = 

In conclusion, each non-zero energy eigenstate with k fermions, will have one and only 
one supersymmetric partner with either k — 1 or /c + 1 fermions. In a space with no holes the 
only eigenstate with zero energy has zero fermions and is the Gibbs measure (See Appendix 
The spectrum is organized as in Fig. 



III. MORSE THEORY 

In the low temperature limit, the organization of the spectrum of allows to derive 
relations concerning the critical points (saddles) of the energy surface, defined as those points 
for which 

N 

|VE|2 = = 0. (24) 

i=l 

Let us study the semidassical low-T spectrum. At the lowest order in T, the potential 
appearing in the Schrodinger operator (jH) 

^='t^\eI (25) 
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is very large except at the critical points of E; so the eigenstates of with low-lying 
eigenvalues are concentrated around the those critical points. 

Let us assume that the critical points are isolated and the Hessian E^ij has non zero 
eigenvalues. Then, as usual, the semiclassical development starts with a harmonic approxi- 
mation around each minimum of W. Consider one of these minima, where the Hessian has 
eigenvalues Ai, . . . , A^. We can develop E in the local coordinates, and going to the basis 
in which the Hessian E^ij is diagonal, we have that, locally: 

Eix') r^Eo + ^xf. (26) 

We can develop H at the first order in T as 

-|A, + TA,afa:|. (27) 

We recognize the Hamiltonian of independent oscillators plus independent fermion 
terms. Along each direction i on the right hand side, we have a harmonic oscillator with 
positive frequency l^jl, plus terms which give —Ai/2 if there is no fermion and +Ai/2 if there 
is a fermion along the direction i. Hence, each fermion term will exactly cancel the zero- 
point energy of each oscillator, provided we have zero fermions if Ai > and one fermion 
if y4j < 0. All in all, we see that we get zero to this order if and only if we have exactly as 
many fermions as unstable directions in the particular critical point 
The next higher eigenvalues A are given by, to leading order: 




A 



i J2 m + 1)1 -A + 2A,n,)] , (28) 



with Ni = 0, 1, 2, . . . and rii is the number of fermions {rii = 0, 1) for each direction. The 
spatial part of the eigenstates are Gaussians times polynomials, thus having widths of order 
^ in the z*'* direction so the approximation is consistent at low T. 
Let us call saddle of index p a critical point whose Hessian has p negative eigenvalues, 
and Mp the number of these. For example the minima are saddles of order while the 
saddles of order A^ are maxima. From the above considerations it follows that around each 
saddle of index p there there is one and only one state with zero energy as T — > 0, and this 
state has p fermions. This means that the Hamiltonian (0)) has, to this order, Mp p-fermion 



eigenstates with a zero energy. As a consequence, 
fermion sector, and the spectrum looks like Fig. |21 



ihere is a gap in the eigenvalues of each 

2a. 
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12 3 4 

l+(3) (3+4) H+2) (2+1) (1) 



FIG. 2: Morse Theory. The arrows indicate the action of Q. The gap in the spectrum means 
that the ratio of the lowest eigenvalue above the gap to the highest below the gap becomes infinite 
at small T. The numbers between brackets indicate the number of states below the gap of each 
fermion number decomposed as in H29|l: the Morse inequalities are evident from the picture. 

Recalling that any non-zero energy eigenstate with p fermions has a degenerate partner 
with either p — 1 or p + 1 fermions (cfr. Section I), one can read from Fig. |21the relations: 

Mo = 1 + 

Ml = i^i + K2, 

Mn = Kn- (29) 

The positivity of the K^s ( -ft'i > Vi) are the strong Morse inequalities. 

We have used the condition that E is defined in a space without holes. As we have seen 
(see Appendix this implies that there is only one eigenstate of zero energy and it has zero 
fermions. If the space has a more complicated topology, there will be several zero-energy 
eigenstates not paired by the supersymmetric charges, and the Morse inequalities become 
slightly more complicated (see Appendix |B|) . 
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FIG. 3: The potential and different eigenstates along the reaction coordinate; (a) is the equilibrium 
density {ijj^^), (b) the first eigenstate - the most stable (tp^^), (c) is the current density in the first 
eigenstate and (d) the spectrum with the gap and the two fermionic sectors. 

IV. STATES AND TRANSITION CURRENTS 

A. A Simple Case 

In the previous section we have given the form of the spectrum, at least to leading order, 
and the corresponding wavefunctions of every fermion number - the latter in the Hermitian 
basis. It may seem that going to the original basis is trivial since it is simply a matter of 
multiplying those approximate wavefunctions by e^^^"^. As mentioned above, this is rather 
tricky, since the factor e^^^'^ will resurrect large deviations which we have neglected. 

Let us first study the simple case of a double well at low temperatures. We consider 
a probability distribution P evolving under the action of the Fokker-Plank Hamiltonian 
(j21) corresponding to an energy as in Fig. El This distribution can be decomposed on 
the eigenstates of the Fokker-Planck Hamiltonian (ie. the zero-fermion eigenstates of the 
Hamiltonian ()13p) 

P(x, t) = Y^ c„^/^°^e-^°*, (30) 

where ip^^ = e~^^ while Aq = 0. 

In the low temperature limit there are two metastable states each concentrated around one 
of the minima. Barrier penetration leads to the Gibbs measure, the symmetric combination 
of those states. In fact the spectrum of the Fokker-Planck Hamiltonian will contain one 
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zero eigenvalue Aq = (the Gibbs measure), one small Ai ~ 0{e~^^^) eigenvalue and the 
rest of them much larger (0(1)). The two pure states, localized on the right and on the left 
are oc ip^^{x) ± ip^^'{x), respectively. If we are interested in the dynamics of the passage 
between the two wells we have to consider times such that the fast relaxation within each 
well has already taken place. At such times, larger than ti ~ j^log{^), we are left only with 
a distribution 

P(x, t > ti) ~ Co^°-^ + ci-^^^e-^'', (31) 

i.e. a combination of states localized to the right and to the left, dependent upon the initial 
condition and time. 

The current at any time is given by: 

^■(-*' - ^ f ) ^<-*) - 5:--^"' ^ §) (32) 

\ / a=l ^ ' 

and its divergence reads: 

(33) 

We can split the contribution of each term as: 

r = A^c^e"^" V"^- (34) 
All the contributions except 1^ eventually vanish, for example: 



at times such that t ^ ta ~ J^^'^di'c 

Equation (jHHj) now implies that the late-time regime current J^{x) is: 



j\x) = ciAie"^^* / dxip^^ix) = cij\x)e-^'\ (36) 

J X 

and this in encoded in the one-fermion right eigenstate 'below the gap' as: 

1^1^) = Q|^i«) = -ijdx j\x) a^\x) ® (37) 

In Fig. 01 we summarize the situation. Note that although the state 1^"'^'*) is sitting on 
the saddle (as we have seen in the previous section), its form in the original basis l^"'"'^) = 
g-/3£;/2|^ih,^^ which encodes the current, is essentially a constant between the two wells. 
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B. States 

We have seen in the previous example how in the low temperature limit one can un- 
ambiguously define metastable states using the eigenstates of the Fokker-Planck operator. 
When the energy function is rough, with several non-equivalent minima, or when the origin 
of metastability is not the low temperature, the construction is less obvious. Suppose the 
Fokker-Planck spectrum has K eigenstates with low eigenvalues 0, Ai, . . . , A^^.i, separated 
by a gap from all the higher ones. One can show that, to the extent that the gap is 
large {\k — ^k-i ^ Ax-i), one can construct exactly K distributions Pq{x), Pk-i{x) by 
linear combinations of the right eigenstates of the Fokker-Plank Hamiltonian ip'^^^x) with 
0<a < K: 

K-l 

P^{x) = 5^T,^^/>^«(^), (38) 

7=0 

such that the Pa{x) are either positive or negligible, and mutually disjoint (the product of 
any two is everywhere negligible), these are the states. We shall take the Pa{x) normalized 
J d^x Pa{x) = 1 V a. Every combination of the right eigenstates 'below the gap' can be 
expressed as a linear combination of the states. In particular, the Gibbs measure is: 

K-l 

r''{x) = J2T,~^'P^{x). (39) 

A useful formula is obtained integrating (jHH|l with respect to x, and noticing that 
(^^oL^^c^R^ = V a > 0: 

T,o = l Va. (40) 

The left eigenstates 'below the gap' are also interesting. By linear combinations of the 
left eigenstates ijj^^{x) with Q < a < K: 

K-l 

A^{x) = J2Ta^i^^\x), (41) 



one obtains functions Aa{x), Ak-i{x) such that each Aa{x) is essentially constant where 
Pa is non-negligible, and is negligible elsewhere. To summarize, right eigenstates below the 
gap are locally Gibbsean, while the corresponding left ones are essentially constant within a 
'state'. The case of very low temperatures is the the simplest one to visualize: in this case 
the Pa are Gaussians sitting each one at the bottom of a local minimum, and the Ai{x) are 
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constant within the corresponding (zero temperature) basin of attraction of each minimum, 
and zero elsewhere. 

The low-lying eigenvalues can also be interpreted as exit times. Under the assumption of 
well separated eigenvalues, their inverses give the exit times corresponding to the metastable 
states. In fact, the low-lying eigenvalues, together with their corresponding eigenstates will 
completely define the long-time dynamics of the system. Indeed, one can define, from a 
probability density P{x), the site-populations as 

= y d^x A^{x)P{x), (42) 

and obtain a master equation for c^{t) as 

dc 

= X] ^^"^^ '^"f" ^ {A^\Hfp\P,,). (43) 

u=Q 

This master equation is accurate at times for which A/<e~'*'^* <^ Ai^_ie~'*'^'"^*. 

Let us conclude by mentioning that there are in cases in which there are more than two 
separated timescales one can make this construction at more than one level - thus obtaining 
different sets of states relevant for the different timescales. 



C. Transition currents from the 1-fermion eigenstates 

Let us assume that the Fokker-Planck spectrum has exactly K eigenstates 'below the 
gap', implying that there are K metastable states. We shall show that, whatever the origin 
of the gap, exactly as in section IIV Al the K — 1 one-fermion partners of these eigenstates 
are the reaction current distributions at long times (for which Xxe.^^^^ ^ Aa'-ic"'^^"^*). 
This may seem rather surprising: the construction of states in the previous section is not 
a priori good in the regions where the probability is negligible. This is not important at 
the level of the probabilities, as those regions carry vanishingly small weight. However, it 
seems to pose a problem for the current, which is important on the barrier - in which region 
the probability corresponding to a state small - and drops to zero within a state, precisely 
where the approximation described in the previous section is reliable. We shall argue now 
that, in spite of this apparent limitations, as we have already seen in the section HV A| the 
1-fermion partners give exactly the long-time currents between states. 
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Let us make the following gedanken-experiment. At time zero we prepare a probability 
density such that it will fall into one state, say tto, < tto < if — 1 . The initial probability 
density can then be written as 

P{X, 0)=Y, Cc.^"'' = cPaA^) + J2 ('^'^) 

a a>K 

We shall study P{x, t) at a time t in the middle of the gap, that is Xxt ^ 1 ^ ^K-it, a 
time sufficiently large so that all the fast components become very small but not yet large 
enough as to populate other states. The density is then 

P(cc,t) = cP„„(a;) + 0(e-^^*) +0(1-6-^^-1*)- (45) 
Let us study now the current 

(46, 

where we have discriminated the contributions to the current of each eigenstate 

(using the notations from the previous section and (P^j) one can see that Cq = cTa^a)- The 
j°'{x) are not normalized, so it is difficult to compare them. In order to do so, we compute 
the divergence of the corresponding terms: 

div r{x) = c„e-^"*A„^/^"^(a;). (47) 

The relative contribution to the current of two terms is of the order: 

where Vq and Va are the regions in which tp""^ > 0, tp^^ > respectively 3|. 

Hence, for large enough times t{\a — \k~i) ^ ln( ^^'^"^^ ^ ), all states a above the gap 
do not contribute to the current. 

In conclusion, we have shown that the escape current of any metastable state is a linear 
combination of the currents associated to states below the gap. This in turn means that, 
within this late-time regime, the current is encoded as a linear combination of some of the 
one-fermion eigenstates 'below the gap', those that have zero-fermion partners 

/K-l K~l 
d^'x Jk{x)al\x) ® I-) = i ^ c^e-^-^ir"") = i Cae'^'^^gi V""")- (49) 

In the next section we interpret those having a two-fermion partner. 



17 



1. Transition times 



Suppose one has the current J escaping a metastable state P(x): 



Mx) oc (T— + eA P{x) ; P{x) = J2 Car'^ix). (50) 



We wish to give an expression for the transition time in terms of the unnormahzed current 
J. For this, we first compute: 



dxi ' / J 1 V dxi 



d^x P e^^ HppP = ^"C^' (51) 



a=l 



and similarly: 



/ d^x e^^^ (div J)2 = I d^'x {HfpP) e^^ {HppP) = X'A 

a=l 



(52) 



because the sums are dominated by the largest eigenvalues Xmax within the sum (j52jl that 
contribute to the state P{x), we have that the smallest escape time is: 

tactiv - \max ' J ^N, ^^E (^iv J)2 " ^^^^ 

Note that the normalization of the current is irrelevant. This formula is valid on the assump- 
tion of separation of timescales, irrespective of its cause. The Kramers expression for the 
low-temperature case can be easily read of this formula, since the numerator is dominated 
by the exponential of the barrier height and the numerator by the exponential of the energy 
of the starting well. If the current is divergence-less (a loop, as we shall encounter later), 
the timescale is infinite. 

An immediate conclusion one draws from ()53|) is that if one knows the current with an 
error 6J{x), it is in the regions with large energy (the saddles) and with large divergence 
(the starting region) where this error translates into a larger error in the timescale. 



D. Loops: blind saddles and subdominant paths 

The one-fermion sector contains in general two kinds of eigenstates states 'below the 
gap'. These are those given by Q acting on a zero-fermion state, and those given by Q 
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FIG. 4: A landscape with a minimum, a maximum and a blind saddle. Below: the low eigenvalue 
spectrum for zero, one and two fermions. The dotted line is the zero level, other eigenvalues are 
exponentially small in l/T. Next higher eigenvalues start at 0(1) (not shown). 



acting on a two-fermion state. The former give us the dominant reaction currents, as we 
have seen already. We now show that the latter give us current loops, and in particular the 
the alternative (subdominant) routes between states. 

The states we are now considering are constructed as follows: given a two-fermion eigen- 
state \p^) = ^l^llPij) ^ ^6 obtain a one-fermion eigenstate as: 



(54) 



unless Ixf ) = 0. From Q\x^) = we immediately conclude that field of the right eigenstate 
is divergenceless: 

(55) 



dxi 



so that if \x ) encodes a single current line, it must be a closed loop. 



1. Simple examples 

A low temperature example will make things clearer. Consider the tilted Mexican hat in 
two dimensions Fig. (j3)): it has a minimum, a maximum, and a 'blind' saddle, one that does 
not lead anywhere. The two-fermion lowest eigenstate is of the form 

\p'') = alal\p^)®\-), (56) 
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where p^{x,y) satisfies: 

- [4. - f ) ^) - (-1 - f ) ^) - ^^"(- (^^) 

which is easily obtained permuting (fermion) particles and holes in the Hamiltonian (fT^ . 
The lowest-lying one-fermion eigenstate is obtained by noticing that the eigenvalue equation 
(j37j) corresponds to the equation satisfied by the left eigenstate of a Fokker-Planck equation 
in a the reversed potential —E{x, y) (cfr. Eq. (f^ ). From the discussion in IIV[ we conclude 
that (the only A{x,y) for the reversed problem) is essentially constant within the region 
spanned by all gradient lines descending from the local maximum (the unstable manifold of 
the maximum, or the stable manifold of the minimum of ~E) - and drops sharply to zero 
at the border of this region. Acting with Q on p^, we obtain the current: 

(x?(-,!/),X?(-,!/))~(^,-|^). (58) 

which is then non-negligible on the gradient paths joining the minimum with the saddle, 
because this is where p^ has a non- negligible gradient. The direction is turnaround, and 
clearly the fiow so obtained is divergence-free. 

Let us now see the general relation between passages and loops with another slightly 
more complicated low temperature example. Consider a situation as in Fig. El There are 
four minima, multiply connected by seven paths going through as many saddles. At low 
temperature, only three of them (shown in thicker lines) have a much shorter passage time, 
and hence dominate the reactions. The other four can be obtained from combination of these 
and the four independent loops — for example one can take each loop going around each 
of the four maxima. The eigenstate structure 'below the gap' refiects this: there are four 
zero-fermion (right) eigenstates corresponding to four minima. One of them is the Gibbs 
measure, the other three have one-fermion partners yielding the three dominant passages. 
The remaining four one-fermion (right) eigenstates correspond to the loops, and they have 
two-fermion partners corresponding to the regions they encircle, including each a maximum. 

The simple tilted Mexican hat problem, and in general any two-dimensional situation, 
allows us also to understand the different roles played by partner eigenstates below the gap in 
the one and two fermion subspaces. As we have mentioned above, all the right two-fermion 
eigenstates which are partners to the loops can be obtained (always in two dimensions) from 
the zero-fermion left eigenstate of the inverted potential. This means that each corresponds 
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1 2 



FIG. 5: A sketch of an energy surface with four minima (full circles), four maxima (open circles) 
and seven pathways passing through one saddles. The thick paths have a low activation times. On 
the right the corresponding spectrum of the Hamiltonian 

to a constant in the region spanned by all trajectories descending from a saddle of index two 
(its unstable manifold), and this will be also true in more dimensions. 



2. Loops: physical meaning and derivations 

The loops have also a physical meaning which may be extended to apply to the nonzero 
temperature situation. Consider a system in equilibrium to which we add a force field hf{x), 
h small, that has only rotational in a restricted region Vr of phase space: 

|^-§^=0 yj,k if xiVr. (59) 

The effect of such a field will be to create currents which will persist even in the stationary 
state. In a system with metastability, these currents can be of two types: those generated 
essentially within a state, and those due to forced passages through barriers; the latter are 
the loops. We shall see that the currents within a state are given by eigenstates above, and 
the loops by eigenstates below the gap of the one-fermion spectrum. 

To make this clear let us go back to the tilted Mexican hat (Fig. H}. Let us consider a 
force whose rotational is concentrated in a restricted 'vorticity' region Vr (the dark region in 
the figure). If the vorticity is concentrated close to the minimum (j3K) the currents generated 
will be due to particles which in a rare excursion happen to fall upon Dp, and then typically 
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FIG. 6: The tilted Mexican hat seen from above. The fuh Unes are the level lines at the energy 
of the saddle (S): the light gray region is above the saddle level. The broken line corresponds to 
the gradient line going from the minimum (m), through the saddle, and back to the minimum 
encircling the maximum (M). The dark gray region corresponds to the domain Pr where the drift 
has a non-zero vorticity. Left: Pr below the saddle level - currents above the gap; Right: Dr above 
the saddle level - currents below the gap. See text. 

fall right back to the state. If we shift the Vr further away from the state, we get a behavior 
of the same kind until we reach a point in which the vorticity is located higher than the 
saddle point, and it becomes more probable for the current to go round the saddle ^p) 
through gradient lines: this is the loop distribution and its essentially independent of the 
exact position of V-p. It is given by the (only) one-fermion eigenstate below the gap. As we 
shall see below, the condition that the vorticity generates a loop around a saddle is that it 
pierces the surface on which the two-fermion eigenstate 'below the gap' is non-zero: this is 
a general fact. 

In order to see this quantitatively, let us study the perturbed Fokker-Planck equation: 

«^. = -E|-(^|- + ^. + ^/.). m 

i ^ ' 

with h a small parameter and /j as above. Proposing a stationary distribution of the form: 

HipPst = ; Pst = Coe-f''' + hP\ (61) 



(co the normalization of the Gibbs measure) the current in the stationary state is obtained 
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as: 



St 



+ E, + hf, ]Pst = h 



OXi 



d 



(62) 



and is zero only if the square bracket vanishes, i.e. if / derives from a potentiaL The current 
is obviously divergence-free, and it is easy to see that the corresponding one-fermion state 
1^^^) = Jd^ 'Ei-Jtix)al\x) ® I-) satisfies: 



0, 



2,hTe-P^\T), 



where 



|r) 



-1 I ^ - — 1 a\al\-). 



dxk dx. 



(63) 
(64) 

(65) 



Equation ()63|) implies that l^"'^) can be developed in terms of one-fermion eigenstates 
with two-fermion partners. Multiplying the Eq. (jMj) by Q, we obtain: 



1^^) = CohH'''Qe-f''\T) = coh H'-'e-^^Q^\T), 



(66) 



where H' is H restricted to the one-fermion subspace. Developing ()66|) in a basis, we find: 

i^"") = c^hj2^\r''){me'')- (67) 

Aq, 

a 

Defining as if+ and the projections of H' above and below the gap, respectively, we 
introduce the components of the current distribution: 



'state 



) = CohH-'e-^^Q^T), 



We can bound: 



(68) 



zohf (r|ge-2^^Qt|r), (eg) 



where A+ mm is the smallest eigenvalue of H in the two-fermion subspace above the gap 
and we have used {ip\AA^ip) < lamaxl"^ ii'li') , with |amax| the maximal eigenvalue of A. As 
(a;|(5^|r) is nonzero only in the rotational region Dp where it can be taken of order one, in the 
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low temperature limit {T\Qe~'^'^^Q^\Ty^'^ ~ e''^^^ where Er = min{E{x)\x G I^r}- Given 
that X^^j^in is of order one and Co is the normalization of the Gaussian, the contribution 
above the gap will bounded as || l^^ate) II— Coe~l^^^ ~ ^-(Er-Emin) ^ exactly what we expect 
of a process that starts in a minimum and climbs up to the region of energy where 
the vorticity is important, and falls back again. 

The contribution of eigenstates below the gap is instead given by the loops, a fact we 
shall show in general in the low temperature limit. The physical meaning of the two-fermion 
wavefunctions below the gap \x°') = Ql^'^^)^ is now clear from equation ()fj7j) : the factor 
(rlQI^""'*) will be important only if the 'vorticity' F intersects the region where the two- 
fermion eigenstate = 01^"''^) is non- negligible. Hence, each two-fermion partner of a 
one-fermion eigenstate below the gap defines the region where a vorticity has to be applied 
in order to excite a current through the corresponding loop. In the low temperature case, 
one expects the current through the saddle to be of order e^^^''°-'^'^^''~^"^^"\ If Ey > Esaddie the 
contribution of the loops dominate, and the bound above means that it can only be given 
by the eigenstates below the gap. Since in simple systems there are only a few of those, the 
distribution will not change dramatically with small changes of the vorticity location. 

In the general case of systems with a gap in timescale, but with nonzero temperature, 
one can still consider forces whose vorticity is 'near' or 'far' from a state, and computing 
the currents induced one can take the construction as a basis for a definition of 'loop'. An 
interesting question is to analyze the effect these loops have in the series development of the 
free energy, an analysis a la Langer would clarify the issue. 

E. Induced currents and holes 

One of the cases in which it is interesting to calculate currents is when we apply a 
constant electric field and join the ends of the sample. We create thus a manifold with a 
hole inside. Up to now we have excluded such situations, and indeed some modifications 
to the arguments have to be taken into account. Consider the simple example of a particle 
in a one- dimensional ring with coordinate x, with < a; < 27r, and no potential. Both the 
Fokker-Planck and supersymmetric operators read: 

H=kQ + Qr = -T^ = HFP, (70) 
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with 

Q = -iT^a ■ Q = -iT^a). (71) 

ox ox 

The zero fermion states are oc e*'^^, with k any integer. In particular, the zero-fermion 
ground state is a constant, ip { expected. On the other hand, one-fermion states 

are also of the form oc e^^^a^\—): we find that we have a one-fermion eigenstate with zero 
eigenvalue \$,^^) = ^o.^ — ), a possibility that we had excluded for spaces without holes (see 
Appendix El)- Furthermore, the one-fermion ground state has no partner: l^"''^) 7^ Q\ij)^). 
The situation changes when we add a constant field so that now: 

H = ^{Q + Qsf = -^(T^ +£) = Hfp, (72) 

with 

Q = -iT—a ■ Qs = -i{T— + £)a). (73) 

The eigenstates do not change, and we still have the same one and zero-fermion ground 
states, but now, remarkably: 

-ir^|l°«)=g^|V°^)^o, (74) 

so that one and zero fermion ground states have become partners. We also conclude that 
the meaning of the one-fermion ground state is to give the stationary current distribution 
around the ring. A last point to see in this simple example is that when the field is on, the 
force does not derive globally from a potential {E = Sx would be multiply valued), and we 
cannot change globally to the Hermitian basis! 

Consider in general diffusion in a space with a hole, so that we can have a force field / 
with everywhere = ^ but not deriving from a global potential. We consider a small 
perturbation 

= - E I- (^1- + ^. + "/.). P«) 

with h a small parameter. Proposing as before a stationary distribution of the form: Pgt = 
Coe~^^ + hP^ the current in the stationary state is again given by (jU^ . It is easy to see that 
the corresponding one-fermion state l^''^) = f d^ic '^iJf {x)a\\x) ®\ — ) now satisfies: 

Qll'') = 0, g||^) = 0. (76) 
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Also: 

\$'') = Qf\Pst) 



\Pst) = -ihe-^''Y.^l 



Cofi + T 



dxi 



(77) 



Just as in the previous example, we have shown that it is the nonconservative field that 
makes IP,*) and l^''^) become partners (otherwise the last of (j77p is empty), and the physical 
interpretation is that each zero eigenvalue one-fermion eigenstate corresponds to a current 
induced around a hole by small fields. Again, for the perturbed Hamiltonian Hpp one cannot 
construct a global Hermitian basis as for Hpp. 

Let us conclude this section with an alternative variational interpretation for the loops. 
Suppose we ask which is the field / such that it maximizes the power W = J /« J/* 

done on the system, while having the minimal Gibbs expectation V for its violation of 
detailed balance (the 'vorticity): 

J d^ic e 



V 

A simple calculation using yields: 



-I3E (dfj_ _ dfk\2 



(78) 



(79) 



where = J d^K fk^k^ and 1^""'^) are 1-fermion eigenstates annihilated by Q (the 
'loops'). In conclusion JF is minimized if f is a left, 1-fermion eigenstate 'below the gap'. 
Clearly, the definition is valid at arbitrary temperatures. 



V. THE BIG PICTURE 

A. Low temperature structures 

At low temperatures eigenstates peak on structures with dimensions smaller than N, 
and fall of exponentially away from them, in a width that vanishes with T. Right eigen- 
states 'below the gap' are made of linear combinations of functions peaked on the following 
structures: 

• Zero fermions: points, the local minima. 
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• One fermion: gradient lines joining minima through saddle points. The eigenstates 
with a zero-fermion partner are the true, open paths, and those with two-fermion 
partners are closed loops. 

• Two fermions: Two-dimensional surfaces containing a saddle with two unstable direc- 
tions, spanned by all the descending gradient lines emanating from it. The eigenstates 
with a one- fermion partner are peaked on open surfaces (surfaces with borders), and 
those with three- fermion partners are peaked on closed surfaces (borderless surfaces). 

The physical property of the open surfaces is that a weak nonconservative field / will 
generate a current turning around their border, and this only if the 'vorticity' region 
in which ( ^ — ^ ) 7^ intersects them. 



• k fermions: /c-dimensional surfaces containing a saddle of index k, spanned by all 
gradient lines descending from it (the unstable manifold of the saddle). 

Again, the eigenstates with k — 1 fermion partners are peaked on open, and those with 
k + 1-fermion partners on closed surfaces. The border of the surface associated with 
the former is the region where the k — 1 fermion partners are peaked. 

Left eigenstates 'below the gap' can be obtained using the fact that a left eigenstate with 
k fermions is a right N — k fermion eigenstate of the problem with the inverted potential 
-E (cfr. Eq. p}))). 

One thus concludes that k fermion left eigenstates are made of linear combinations of 
functions peaked on the following structures: 

• Zero fermions: constant within a basin of attraction of each minimum (a well known 
fact). 

• One fermion: N — 1 dimensional basin of attractions of saddles (themselves subsets 
of the borders between basins) . 

• Two fermions: N — 2 dimensional basins of attraction of saddles with two unstable 
directions. 

• k fermions: N — fc-dimensional surfaces spanned by the set of descending paths ter- 
minating in each saddle of index k (i.e., the basins of attractions of these saddles). 
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This is the structure of basins within basins that was argued is relevant in systems with 
slow dynamics 3|. One can again distinguish open and closed surfaces, and this is related 
to whether the wavefunction of k fermions has a partner with k — 1 ot k + 1 fermions. 

Each of the right eigenstates below the gap are peaked on unstable manifolds of the 
corresponding critical point, while the left eigenstates are peaked on the stable manifolds. In 
this frame, the Q and Q operators act as 'boundary operators' generating what is called in the 
mathematical literature the Morse (co)homology. We shall derive these results constructively 
in the next section. 

A simple, three-dimensional example will illustrate this. Consider the energy E{x, y, z) = 
e{x)+e{y)+e{z), where e{xi) is the symmetric double-well potential of Fig. IHl The landscape 
is defined by a cube with minima in its 8 vertexes, 12 saddles of index 1 midway along its 
sides, 6 saddles of index 2 at the face centers and one saddle of index 3 at the cube's center. 
The right eigenstates below the gap are as follows i) 8 with zero fermions located at the 
vertices's, ii) 12 with one fermion located on the sides, of which six are passages and six 
are loops along the perimeter of the faces, in) 6 with two fermions peaked on the faces, of 
which five are independent open faces and one is the total (closed) surface of the cube, and 
iv) one with three fermions constant inside all the interior of the cube. 



B. Defining 'free energy' structures 

The zero-temperature limit allows us to see very clearly the different structures that 
emerge. However, the main point of this paper is that these can be transferred to a more 
general situation, provided that there is timescale separation - whatever its origin. As we 
have seen already, the role of local minima in the low-temperature situation is taken by 
metastable states, and the role of gradient lines by reaction current distributions. In the 
previous section we also attempted a general definition of 'reaction loop', on the basis of the 
currents that can be induced by a non-conservative weak force. One has the possibility of 
also defining higher structures associated with fermion subspaces of higher fermion numbers, 
(like borders, basins etc.) in general. In problems in which one can define a free-energy 
landscape in a precise way (typically mean-fieldish situations), the structures we have defined 
should recover their geometric appearance: states becoming points, reaction distributions 
becoming lines, etc, but now in the free-energy landscape, in which each point stands for 
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many configurations. The importance of tfie construction we liave been describing is tfiat 
it does not rely on sucli a landscape: the only assumption is timescale separation. Once 
this is given, the Morse-theory constraints on the objects follow automatically, thus showing 
that the construction makes geometric sense. We shall also see in what follows how these 
structures are be approached by higher forms of stochastic equations. 



In order to obtain the currents we must find eigenstates of the supersymmetric Hamilto- 
nian (fT^ 'below the gap' in the 1-fermion sector. This may be achieved by solving © at 
times larger than the microscopic times starting from several initial configurations. In the 
zero-fermion case, one in fact simulates the Langevin dynamics (0) rather than solving the 
Fokker-Plank equation, in order to obtain metastable states. The question naturally arises 
as to which is the diffusion equation that reproduces (jH)) in the one fermion subspace, and 
more generally in any /T-fermion subspace (3fll |. 

Let us do this for one- fermion wavefunctions first. 

Consider first one particle carrying an A^-component vector degree of freedom u. Let the 
position of the particle evolve as a Langevin process (^, and the vector u as: 



VI. DIFFUSION DYNAMICS 



A. Dynamics 




(80) 



3 



From and (IHUj). we have that the joint distribution function F{x,u,t) evolves then as: 




(81) 



Consider now the evolution of the partial averages: 




(82) 



dRa{x,t) 
dt 




'a 



dF{x, u, t) 
di 




a 



Hfp + ^(*' ^' ^) 



'i 



HppRa{x, t)-J2 EajRjix, t), (83) 



j 
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i^IG. 7: Snapshot of a population of walkers in the stationary state. The potential is taken from 
3^ . It has two minima (right and left), two saddles (top and bottom) and a maximum in the 
center. 

where we have integrated by parts on Mq. This is equation (0), as announced. 

The evolution as in (j8U|) and (j81|) has the problem that the particle may very rarely 
visit a given region of space, but once it does so have large components for u. A practical 
modification is to preserve the norm of the vector attached to the particle. Putting v = u/\u\ 
we obtain the following equation for a function F{x, v, t): 



dF{x,v,t) 



Hpp + ^-^-1 -'"i^ VkViEkl j> - ^ VkViE, 



ij (. kl J kl 



dt 

Computing the evolution of the partial averages 



kl 



F{x,v,t). (84) 



Ra{x,t) = I dvVaF{x,v,t), (85) 

the result is again (|7j). The diffusional process involved is however quite different. The first 
term in the square bracket in (jS^ tells us that the dynamics of the particle is still of the 
Langevin form. The second term now gives for the evolution of v: 

Vi = - ^ Eij{x)vj + Vi ^ VkViEki, (86) 
j kl 

which preserves the condition || v \\= 1. The third in (j84j) is a 'cloning' term, creating and 
destroying particles at a rate Ylki'^k^iEki 0]. In Fig. [7|we show the numerical solution of 
the diffusion equation for vector walkers in a potential. 
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B. Low temperature structures 

Let us now use the equations for v to show that at low temperatures the eigenstates with 
eigenvalues close to zero correspond to lines joining minima through the saddles. At low 
temperatures, particles fall along gradient lines: 

X, = -E,i{x). (87) 

Assume that the distribution consists of particles sitting along an isolated gradient line, 
and having Vi = Ei/\VE\. Putting this into ()8fi|l . we find that the condition Vi = Ei/lVEl 
is preserved as particles fall along the gradient line. On the other hand, the depletion 
of particles along the trajectory consists of a term due to migration — ^ = '^^"^'^ {di 
the element length along the line), and of cloning Ylki'^kViEki: both terms cancel since 
di is parallel to v. Hence, a uniform distribution of particles along a gradient line with 
Vi = Ei/\VE\ is stable, provided nothing happens at the ends. Now, the only possibility 
for the ends not to destroy stability is that they are stationary points, so that there is 
no particle exchange there. Furthermore, the distribution has to be in particular peaked 
along a path joining two minima through a saddle of order one. The reason is as follows: 
particles are constantly falling, the measure is preserved because there is a high birth rate 
near the saddle. Now, if the saddle in question is of index higher than one, the slightest 
noise will make particles that are born near it emigrate in other directions, as there is more 
than a single descent path in that case, rendering the solution along a single gradient line 
unstable. A stable one-fermion solution on a higher dimensional surface is on the other hand 
impossible because the surface expansion rate is no longer compensated by the cloning. 

In Appendix O we give the generalization of the evolution equations for higher fermion 
numbers, and we use them to generalize this argument to show the result announced in 
the previous section that low-temperature fc-fermion eigenvalues below the gap are linear 
combinations of constant densities filling the surfaces spanned by all the descending paths 
emanating from saddles of order k. The argument is entirely similar to the one for the one- 
fermion sector: each particle has a /c-form attached to it, whose distribution is preserved 
as it falls down. The cloning term precisely compensates the effect of the redistribution of 
particles, and the solution for k forms is unstable unless particles are born near a saddle of 
index k. 
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C. Long-time evolution of v 

At times longer than the inverse of the lowest one-fermion eigenstate, one can expect that 
F{x^ v,t) will converge to a stationary distribution symmetric in v. 

lim F{x, V, t) = F^'^ix, v), F^'^{x, v) = F'"^{x, -v) \/v, (88) 

t— >oo 

and the averages (jH3j) vanish 

Rl%x, t)= I d% Va F^'^ix, V, t) = 0. (89) 

J|ll|2 = l 

To compensate this mean death rate, the usual practice in Diffusion Monte Carlo schemes 
is to add an overall cloning probability (see js^). In any case, one can show that F'^'^{x, v) 
itself will be peaked paths at low temperatures. To do this, it suffices remake the argument 
of Section El A. 

VII. PATH SAMPLING: LYAPUNOV WEIGHTS 

The same ideas can be written in the path-integral formalism. Let us start by computing 
J(a;o,a;i) = ^(-| ® {xo\ai e~^^ ajl^i) ® I") = 

i 

= E^r'(^i)ef (^o) e-'-\ (90) 

ai ai 

Because H is quadratic in the fermions, the evolution for the a] is linear, and we have in 

nn 

terms of the trajectories |33l. I34| : 

/(xo,xi) = Q D[paths] TiVpath S[xi + Ej - r]i]j = J D[paths] TiVpath e-^'"'*^ 

(91) 

where the sum is over all paths going from x -i to Xq, and the average is over the noise rj 
realization. Spath is the usual Langevin action [23 

U(t) is the matrix solution of the linear equation 

Ui, = -Y,E,kUkj UiO) = I, (93) 
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FIG. 8: One typical trajectory with the ends fixed in two difi^erent metastable states sampled with 
Lyapunov+Langevin action. Same potential as in Fig. [7] 

which depends on the path through Ei^. It describes the linear transformation of a small 
region around the trajectory, defined by a set of nearby initial conditions and the same 
thermal noise. With these notations, we have: 

Iixo,Xi) = J^^fM^fMe-^"' = J D[paths] e-^^^'^'^-^'^-^^\ (94) 

ai 

where we have defined the (pseudo) Lyapunov exponent (with the time included!), for large 
t: 

Ll^,, = In [TrU,.,,] . (95) 

The prefix 'pseudo' is a reminder of the fact that actually, true Lyapunov exponents is 
defined on the basis of the trace of UU^ For Xq = X\ our matrix U is symmetric on average 
{U'ij) = {Uji) (a consequence of detailed balance), but not along a single trajectory. We shall 
return to this point later. For large t, L^p^th = l-^maxl where \^ax is the eigenvalue of \Jpath 
with the largest real part. Note that L^path can be calculated for long times by considering 
the path and a nearby path starting from an initial condition close to Xq and evolving with 
the same noise, just as in the computation of an ordinary Lyapunov exponent, or on the 
basis of the force required to keep the distance between paths fixed. 

Consider the trajectories weighted with the modified action (jHH). We wish to know the 
distribution of x at an intermediate time t' . Let us compute the expectation of an arbitrary 
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FIG. 9: One typical closed trajectory sampled with Lyapunov+Langevin action. Same potential 
as in Fig. [7J 

function B{x) at t': 

{B{t', xo, xi)) = r\xo, xi) I d^x ^(-1 ® {xo\aie-^'-''^" B{x)e-''^ a\\x-i) ® |-) 

= I-\xo,x^) I d^xY,if{x^)^t''{xo) {e''\B{x)\i^'')e-^P^'-''^-^-r (96) 

Considering t' and (t — t') longer than the microscopic times but much shorter than the 
passage times one can retain only the contributions 'below the gap': 



{B(t\xo,x^)) = I-\xo,x^) /"rf^ 5^C^(a.i)CVo) (r^lSNII^^) 



(97) 



x] 



Now, if Xq and x\ are at the states at the ends of a reaction, the factor ^f^(a;i)^f "^(ajo) 
'^.i^f^{xx)^f^{xo)e~^^^^'^'^ selects the relevant currents , and we obtain: 

(S(t',Xo,a;i))~ j d^'x {e\B{x)\i'') = J d'^x {e\B{x)\e) = J d'^x \\ ^i^) f B{ 

(98) 

where we have assumed the reaction is given by a single 1^"'^). What we have shown is that 
long paths sample the barrier || $,'^{x) \\. In other words, trajectories have ends of the order 
of the microscopic time in Xi and Xq, but otherwise spend most of their time in the barrier: 
see Fig. |H1 If we consider closed paths without restrictions on the starting point, we have: 

d%{B{t',Xo,Xo))= I d%d''x Y.if{Xo)it''{^o){C''\B{x)\e'') 



af3i 



^ I d^x {C^\^'^^){C^\B{x)\i'^^) = J2 f d^x \^f{x)\^B{x) 

af3 ia 



(99) 
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FIG. 10: One typical trajectory with one end fixed in a metastable state and the other left free, 
sampled with Lyapunov+Langevin action. Same potential as in Fig. [Tj 

where the sums are only over the 'eigenvalues below the gap'; this means that we perform 
a fiat sampling over all barriers: see Fig. IHl If instead of || ^'^(jc) ||^ we wish to sample the 
squared current Yli^i^^i^ ~ II $^{^) 11^? we have to add a Gibbs weight at a single time 
in the closed path measure. 

We can of course fix only one end, and the situation obtained is as in Fig. El sampling 
the escape paths from one metastable state. 




A. Higher index barriers 



The procedure outlined above can be generalized in a straightforward way to higher 
indices. One starts from: 



(100) 



which will select the fc-fermion eigenstates 'below the gap'- Again, the evolution of each 
fermion is linear, and a straightforward calculation j^J, shows that the path-integral 
reads: 



(101) 



where i^^/j is the (pseudo) Lyapunov exponent, defined as 



''path 



In 



detp{U, i\ 



(102) 



Ui,---,«fe 
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where detp{U, ii, . . . , ik) are the k minors of U. 

""path 



For large t, L^^^^t^ is the logarithm of the sum of the k eigenvalues Upath having the largest 



real parts. L^^^^^ is a measure of the expansion of /c-dimensional surfaces defined by nearby 
trajectories subjected to the same noise Q]. Just as before, trajectories weighted with L'^p}^y^ 
will pile up in index k barriers. 

Before concluding this section, let us point out that, unlike the true Lyapunov exponents 
defined in terms of U^U, the ones we are using here are real on average, but there could be 
rare trajectories for which they are imaginary. The practical procedure is then to separate 
trajectory space in those that have a real, and those that have imaginary value of Lp^^^^. As 

(k) 

usual in these cases, the separation is natural since LpJ^/^ diverges in the frontier. 

VIII. CONCLUSIONS 

In this article we have shown how the constructions based on supersymmetry can shed new 
light on statistical mechanical questions, providing definitions and computational schemes 
for barriers beyond the low-temperature or the mean-field cases. 

We have deliberately avoided maximal generality at each step, in various cases leaving the 
most general derivations for the Appendices. We have also not attempted full mathematical 
rigor. The aim has been to convince the reader that all developments are elementary, though 
we believe quite useful. There are three important omissions: 

• Continuous symmetries, leading to non-isolated saddle points and barriers: the subject 
of degenerate Morse theory 17[. The formalism adapts itself rather naturally to this 
case, so we are confident that the discussion in this direction can be made more 
complete. 

• Dynamics with inertia (Kramers equation). This is important for practical applica- 
tions, in which reaction paths have to be found in systems with inertia. 

Hamilton's equations do possess a supersymmetry, as shown by Gozzi and Renter 
(see |38| and references therein), who used it to rederive some very early results by 
Ruelle where Hamiltonian, as opposed to Langevin, dynamics is used to study the 
topology of the space. Part of this supersymmetry survives for the Kramers equation 
40|, so the results in this paper, and indeed all the construction related to Morse 
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Theory, can be extended to that case. This is a promising hne of research, with a 
considerable number of apphcations including, apart from the study of reaction paths, 
the determination of Ruelle-PoUicott resonances in chaotic systems. 

One can ask if the approach presented in this paper is applicable only to smooth 
energy functions in space, since there are real problems with energy functions which 
have singularities: exclusion or Coulomb interactions, hard walls, etc. It is possible 
to obtain rules for the behavior of particles (or rather, the forms attached to them) 
after a collision with a hard wall without having to integrate the bounce trajectory 
every time, by deriving the effect of a regularized wall in the limit of infinite steepness. 
Using this method it is possible to construct diffusion equations for problems which 
are entirely entropic, such as hard spheres. 

The question of the application to full quantum evolution remains open. There is of 
course the less general possibility of coupling the present scheme to a Carr-Parrinello 
approach, the dynamics being essentially classical in that case. 







Work is in progress (40, |4]|, stimulated by the prejudice that things that are pleasant 
should also be useful. 



Acknowledgments. We thank A. Oancea for useful discussions. 

APPENDIX A 

In his Appendix we sketch the proof that the only zero energy eigenstates of the Hamil- 
tonian ()13|) is the Gibbs measure (PT|). The proof will be made using the hermitian basis. 

Each eigenstate of energy zero must be annihilated by both Q and Q. Indeed, the 
definition of the zero-energy state: 

+ g'^^g'^)!^^'') = ^(11 Q'^m f + w Qm in, (m) 

is equivalent to 

g'^IV'^) = 0, g^lt/;'^) = 0. (A2) 
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Let |t/> ) be a state with p > fermions such that Q ) = 0. If we construct the state 



-'0 



one can easily verify that 



This in turn means that 



xiQ'^m = mQv) =11 m r> o, 



(A3) 



(A4) 



(A5) 



which is incompatible with Q'^\ip^) = 0; we have proved that states with 1 fermion or more 
cannot be annihilated simultaneously by and Q'^ therefore they cannot have zero energy 



42|. 



(A6) 



On the other hand, if has no fermions, ()A2|) can be written as: 

and this has only one solution ilj'^{x) = CqC"'^'^^^-'/^ (E^- 

Using standard arguments one can show that, for an energy E bounded from below and 
satisfying (jS)), \x^) ^iid e~^^/'^ have a norm, and thus are in the Hilbert space associated 
with H^. The conclusion is that the only zero energy state of is 



-I3E/2 



(A7) 



APPENDIX B 



In this appendix we write the Morse inequalities derived in section IIIII for the trivial 
topology in a more general context. The number of exact zero energy states, for each 
fermion sector (let us call them Bp) does not depend on the energy but only on the topology 
of the space. To see this, suppose that the energy is changed by -E'(x) -E'(x) + 6E{-x). A 
short computation yields, to first order: 



= 

2T 



1 

2T 



fBll 



and this has zero matrix elements between states with zero eigenvalue, as they are annihilated 
by the charges. First order perturbation theory tells us then that the eigenvalues stay zero. 
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One can also exclude the possibility of a non-zero eigenvalue becoming zero, by applying 
the previous argument to the reverse perturbation. 
For , we have shown that the Bp are 

Bo = l, Bi = 0, Bn = 0. (B2) 

Thus, following the same arguments as in section llTTl one can write generalize the equalities 
^ to: 

Mo = Bo + K,, 

Ml = B^ + + K2, 

Mn = Bn + Kn- (B3) 

where, again, Ki > 0,Wi. 



APPENDIX C 



The evolution equation ^ 



for a vector 



(CI) 



(C2) 



reads, in components: 



{^2),■■■,a,...,a{^^.)^ 



(C3) 



where a denotes all permutations of k indices, and n{a, a) is the sign of the permutation 
{ii,i2, . . . ,a, . . . ,ik) (a, cr(z2), • • • , (^{ik))- The tpa{h),...,a{ik) are antisymmetric with 
respect to permutations of indices. Proposing the evolution for functions of F{x, u, t), where 
u is the set Ui-^^,,,^i^_, themselves completely antisymmetric: 



dF 
Itt 



d 



FP 



(i2),...,a,...,o-(ifc) 



Fix,u,t), (C4) 
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we can check integrating by parts that ipi-^^^^^^i^lx) = J d^i Ui^^,,,^i^ F{x, u, t) evolves according 
to ()(]3j) . This in turn means that x evolves according to the Langevin equation, while 



We can also write an equation for normalized variables: 

Ui 



"It,..., Ik 



such that (IC4I) becomes: 



dF 



(C6) 



+ M{v) 

where: 



ii2),...,a,...,a{ik) - Vh,-,ik ■^i'") 

Fix,v,t), (C7) 



ii,...,ik o" ct 

The particles perform Langevin diffusion, while the equation of motion for the v read: 

Vi,,...,ik = -J2Yl (-1)"^"'"^ E^ihU ^^<xfe),...,a,...,<xfc) - Vi,_i, Afiv), (C9) 

thus preserving the normalization ()C6|) . There is also cloning, proportional to A/'(v). 

The equations of motion for v have an interesting interpretation. Consider a point a; and 
a small oriented /c-volume element determined by 

V'' = {x + 5xi)A{x + Sx2)A---A{x + 5xk) = M ^ ?;i,,.„,,,ei, A ■ ■ ■ A e^,, (CIO) 

where A is the external (wedge) product and ei are the basis vectors. We have separated 
the norm Ai from the (normalized) 'direction' v of V'^. It is straightforward to see js^ 
that equation ()C9|) indeed gives the evolution of v as the points are carried by the drift, and 
M = M. gives the expansion rate of the norm of . This property is at the basis of the 
use of the present formalism to study Lyapunov exponents |34| . 

Now we can outline a proof that the right /c-fermion eigenstates 'below the gap' are 
concentrated on /c-dimensional surfaces spanned by the trajectories descending from a saddle 
of index k. Let us propose that on such surface we have particles whose v at each point is 
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tangential to such a surface (i.e. it can be generated as in ()C10|) ). and the density of such 
particles is constant along the surface. Both features are preserved by the evolution. First, as 
the particles go downhill, their forms change so as to remain tangential: this is because their 
evolution are precisely based on the linearized equation on the tangential space. Secondly, 
the cloning rate is exactly opposite to the expansion rate of a small volume advected downhill: 
as a region expands its population increases and vice-versa. Furthermore, close to the saddle 
point the surface density expansion rate is —Ai, . . . , —A^, {Ai the Hessian's eigenvalue) while 
the cloning rate for a p-form is: —Ai, . . . , —Ap. Because by assumption Ai < 0, . . . , A^ < 
and Ak+i > 0, . . . , A^ > 0, the cloning is insufficient to maintain a stationary situation 
unless p = k. 
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